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Abstract 

Nonlinearity  and  the  time- varying  dynamics  of  fuel  cell  systems  make  it  complex  to  design  a  controller  for  improving  output  performance.  This 
paper  introduces  an  application  of  a  model  reference  adaptive  control  to  a  low-power  proton  exchange  membrane  (PEM)  fuel  cell  system,  which 
consists  of  three  main  components:  a  fuel  cell  stack,  an  air  pump  to  supply  air,  and  a  solenoid  valve  to  adjust  hydrogen  flow.  From  the  system 
perspective,  the  dynamic  model  of  the  PEM  fuel  cell  stack  can  be  expressed  as  a  multivariable  configuration  of  two  inputs,  hydrogen  and  air-flow 
rates,  and  two  outputs,  cell  voltage  and  current.  The  corresponding  transfer  functions  can  be  identified  off-line  to  describe  the  linearized  dynamics 
with  a  finite  order  at  a  certain  operating  point,  and  are  written  in  a  discrete-time  auto-regressive  moving-average  model  for  on-line  estimation  of 
parameters.  This  provides  a  strategy  of  regulating  the  voltage  and  current  of  the  fuel  cell  by  adaptively  adjusting  the  flow  rates  of  air  and  hydrogen. 
Experiments  show  that  the  proposed  adaptive  controller  is  robust  to  the  variation  of  fuel  cell  system  dynamics  and  power  request.  Additionally,  it 
helps  decrease  fuel  consumption  and  relieves  the  DC/DC  converter  in  regulating  the  fluctuating  cell  voltage. 

©  2008  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Proton  exchange  membrane  fuel  cells,  also  known  as  poly¬ 
mer  electrolyte  membrane  fuel  cells,  have  the  advantages  of 
high  power  and  energy  density,  low  operation  temperature,  fast 
start-up,  and  higher  efficiency,  but  still  face  many  challenging 
technical  hurdles,  such  as  cost,  size,  weight,  durability,  reliabil¬ 
ity,  and  air,  thermal,  and  water  management  [1-3].  In  addition, 
fuel  cells  are  subjected  to  various  situations  of  time- varying  load, 
during  which  the  air-flow,  gas  pressure,  temperature,  humid¬ 
ity,  and  membrane  hydration  must  be  controlled  over  a  wide 
range  of  operation.  Typical  fuel  cell  system  operation  requires 
prompt  measurement  of  system  states  by  a  set  of  sensors,  such 
as  flow  meters,  thermocouples,  pressure  transducers,  voltmeter, 
hall  sensor,  hydrogen  detector,  etc.  These  signals  are  typically 
fed  back  to  a  microprocessor  that  calculates  proper  control 
actions  based  on  a  specific  control  strategy  that  is  executed  by 
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various  actuators,  such  as  an  air  pump,  a  humidifier,  solenoid 
valves,  fan  motors  and  safety  devices. 

The  complex,  nonlinear  dynamics  of  fuel  cell  systems  are 
usually  approximately  described  by  the  principles  of  electro¬ 
chemistry,  fluid  dynamics,  and  heat  transfer  in  terms  of  physical 
parameters,  material  properties,  and  universal  constants  under 
various  assumptions  and  constraints.  The  dynamic  response  of 
cell  voltage  and  current  involve  reactant  pressure  and  flow  rate, 
local  species  concentration,  and  load  electrical  demands,  which 
can  vary  rapidly  enough  that  many  assume  these  are  instanta¬ 
neous.  On  the  other  hand,  fuel  cell  temperature  and  its  regulation, 
humidity,  cell  hydration,  as  well  as  overall  heat  management  pro¬ 
cesses  are  characterized  by  relatively  slow  dynamic  response,  as 
analyzed  by  Mueller  et  al.  [4].  Many  researchers  have  devoted 
themselves  to  developing  and  using  steady-state  models  of  PEM 
fuel  cells,  describing  the  relationships  amongst  physical  vari¬ 
ables  through  the  Nernst  equation,  gas  diffusion  equation,  and 
species  concentration  and  polarization  equations.  Most  recently, 
an  increasing  number  of  researchers  have  focused  on  dynamic 
models  to  describe  the  transient  response  of  PEM  fuel  cell  sys¬ 
tems.  Pukrushpan  et  al.  [5]  modeled  transient  dynamics  with 
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a  set  of  first-order  differential  equations  governing  the  air  com¬ 
pressor,  mass  transportation,  energy  conservation  of  the  reactant 
flows,  and  pressures  in  the  cathode  and  anode.  Ceraolo  et  al.  [6] 
provided  more  precise  partial  differential  equations  to  describe 
both  static  and  dynamic  behaviors  of  the  PEM  fuel  cell,  including 
gas  diffusion,  proton  concentration,  and  mass  transportation.  In 
their  spatial,  time-dependent  fuel  cell  model,  Golbert  and  Lewin 
[7]  included  the  dynamics  of  water  condensation,  evaporation, 
and  generation,  as  well  as  quasi- steady- state  temperature  profile. 
Pathapati  et  al.  [8]  derived  a  dynamic  model  with  the  effects  of 
charge  double-layer  capacitance,  dynamics  of  flow  and  pressure, 
and  mass/heat  transient  features  of  fuel  cells,  which  predicted  the 
transient  response  of  cell  voltage,  temperature,  gas  flow  rates, 
and  pressure  under  a  sudden  change  in  load  current.  A  highly 
dynamic  PEM  fuel  cell  stack  model  was  proposed  by  Shan  and 
Choe  [9],  which  took  into  account  the  most  influential  property 
of  temperature  on  the  output  performance  and  dynamics. 

Although  some  researchers  have  proposed  an  effective  model 
for  system  and  control  development  [4,10],  the  nonlinearity 
and  time-varying  characteristics  still  pose  difficult  problems 
for  system  identification  and  control.  Simplified  models,  either 
linear  or  nonlinear,  with  nominal  parameters  are  usually  valid 
within  a  linear  range  of  operation.  Such  models  are  usually  used 
to  investigate  stability,  sensitivity,  observability,  and  control¬ 
lability,  before  designing  a  controller.  More  practically,  these 
parameters  that  vary  with  time  and  operating  states  are  identi¬ 
fied  in  real  time  to  match  system  response  with  minimal  errors. 
Therefore,  adaptive  control  can  serve  as  a  feedback  law  to 
achieve  control  objectives  subject  to  the  variation  of  system 
parameters  as  well  as  external  disturbances.  Pukrushpan  et  al. 
[5]  designed  an  observer-based  feedback  controller  to  protect 
the  fuel  cell  stack  from  oxygen  starvation  during  the  change  of 
current  command,  using  the  linear  quadratic  technique  based  on 
the  linearized  state-space  model.  Paradkar  et  al.  [1 1]  integrated  a 
linearized  PEM  fuel  cell  model  into  a  power  plant,  and  simulated 
a  load  frequency  control  problem  by  an  optimal  controller  based 
on  the  disturbance  accommodation  control  theory.  Schumacher 
et  al.  [12]  employed  a  fuzzy  controller  to  manage  humidity  by 
adjusting  the  fan  voltage  for  supplying  air  to  a  miniature  PEM. 
Jurado  and  Saenz  [13]  noted  that  their  proposed  adaptive  control 
strategy  was  able  to  stabilize  a  hybrid  fuel  cell  and  turbine  system 
subject  to  system  parameter  variations  and  external  disturbances. 
Tekin  et  al.  [14]  provided  a  fuzzy  proportional-derivative  (PD) 
controller  plus  integral  (I)  controller  to  adjust  the  air-flow  rate 
of  a  5  KW  PEM  fuel  cell.  It  was  found  that  the  ratio  between 
the  consumed  energy  by  the  motor-compressor  group  and  the 
energy  delivered  by  the  fuel  cell  was  greatly  reduced.  Woo  and 
Benziger  [15]  applied  a  PID  controller  to  regulate  the  output  cur¬ 
rent  of  a  single-cell  PEM  fuel  cell  by  controlling  the  flow  rate 
of  hydrogen,  while  the  full  utilization  of  hydrogen  was  achieved 
with  30%  stoichiometric  excess  on  the  oxygen. 

Recently,  Meyer  and  Revankar  [16]  surveyed  the  control- 
oriented  models  and  control  strategies  of  PEM  fuel  cell  systems. 
The  researchers  indicated  that  future  attempts  to  develop  nonlin¬ 
ear  multi-input  multi-output  systems  must  be  important  over  the 
entire  operating  range  in  order  to  seek  reduced  reactant  usage 
under  various  power  demands.  Methekar  et  al.  [17]  proposed 


a  linear  multivariable  model  for  controlling  the  average  power 
density  and  solid  temperature  of  a  PEM  fuel  cell  to  avoid  oxygen 
starvation.  The  selected  input  variables  were  the  inlet  molar  flow 
rate  of  hydrogen,  the  inlet  molar  flow  rate  of  oxygen,  the  inlet 
flow  rate  of  coolant,  and  the  inlet  gas  temperatures  of  anodes  and 
cathodes.  These  researchers  simulated  the  order  reduction  from 
a  distributed  parameter  model  to  a  set  of  first-  or  second-order 
transfer  functions,  by  which  the  Pl-controller-based  power  and 
solid  temperature  control  loops  were  constructed. 

The  authors  have  also  presented  a  preliminary  study  of  system 
identification,  and  robust  and  adaptive  controls  of  a  low-power 
PEM  fuel  cell  [18,19].  This  paper  extends  the  previous  single¬ 
input  single-output  (SISO)  adaptive  control  to  a  multivariable 
configuration  of  indirect  model  reference  adaptive  control.  Both 
the  air  and  hydrogen  flow  rates  are  used  to  regulate  cell  voltage 
and  current  according  to  the  updated  information  of  the  external 
load  and  system  dynamics. 

2.  Material  and  methods 

2.7.  Fuel  cell  specifications 

The  fuel  cell  stack  in  use  consists  of  15  cells  with  an  active 
area  of  about  50  cm2,  as  shown  in  Fig.  1 .  The  cells  are  connected 
with  a  pre-treated  membrane  (Nation®  1 1 2)  by  hot  press  for  opti¬ 
mum  conditions,  and  are  electrically  connected  in  series  with 
catalyst  layers  where  the  platinum  coating  is  about  0.2  mg  cm-2 
at  the  anode  and  0.4  mg  cm-2  at  the  cathode.  The  gas  separa¬ 
tor  plates  are  made  of  the  graphite  material  (TOYO  TANSO, 
Japan)  treated  to  become  gas  impermeable;  each  is  as  large 
as  160Lmm  x  50  W mm  x  2  T mm  with  two  fins  on  each  side, 
and  each  fin  has  a  size  of  160Lmm  x  15  W  mm  x  17  mm.  This 
design  is  so  similar  to  the  cooling  fins  of  internal  combustion 
engine  that  it  increases  the  cooling  surface  and  improves  the 
cooling  efficiency.  In  addition,  the  hydrogen  flow  field  consists 
of  a  serpentine  pattern  of  parallel  double-channels,  and  the  air¬ 
flow  field  consists  of  a  straight  pattern  of  parallel  multi-channels 
and  ribs  in  an  active  area  of  40  mm  x  130  mm. 

The  rated  and  peak  power  outputs  of  the  fuel  cell  are  1 17  W 
at  9  V  and  124W  at  7.9  V,  respectively.  The  maximum  fuel- 
to-electricity  efficiency  that  the  fuel  cell  can  achieve  is  about 
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Fig.  1.  100W  fuel  cell  stack. 
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Fig.  2.  PEM  fuel  cell  system  dynamics. 


37%  on  a  lower  heating  value  (LHV)  under  the  dry  H2/air  and 
humidification-free  conditions.  The  fuel  cell  system  was  a  pro¬ 
totype  from  Delta  Electronics.  The  whole  system  includes  an 
individual  temperature  control  loop,  water  management  loop, 
and  flow  control  loop.  This  paper  intends  to  replace  the  flow 
control  loop  using  adaptive  control  while  keeping  the  others 
unchanged.  The  flow  rate  of  natural  air  is  adjusted  by  an  air 
pump,  while  the  flow  rate  of  hydrogen  is  regulated  by  a  solenoid 
valve. 

2.2.  Fuel  cell  system  model 

In  order  to  develop  the  control  strategy  of  the  low-power  fuel 
cell,  the  dynamic  models  from  the  previous  studies  are  integrated 
to  describe  the  mass  and  heat  transfer  in  flow  channels,  Nernst 
equation,  cathode  gas  diffusion,  cathode  kinetics,  and  proton 
concentration  dynamics  [6,20-23].  From  modern  control  theo¬ 
ries  [24] ,  the  control-oriented  fuel  cell  system  block  diagram  can 
be  approximated  as  a  two-input  two-output  system,  as  shown  in 
Fig.  2.  Equipped  with  additional  balance-of-plant  (BOP)  com¬ 
ponents,  such  as  an  air  pump  and  hydrogen  valve,  the  fuel  cell 
system  can  be  simply  reduced  to  a  standard  block  diagram,  as 
shown  in  Fig.  3. 

From  the  system  viewpoint,  the  hydrogen  at  a  certain  pressure 
is  fed  by  a  solenoid  valve  whose  dynamics  are  denoted  by  the 


Fig.  3.  Approximate  fuel  cell  system  block  diagram. 


block  Gs  that  describes  the  relationship  between  the  hydrogen 
flow  rate  Ah  and  valve  voltage  Vs .  On  the  other  hand,  air  at  the 
standard  atmosphere  of  1  bar  is  often  transported  by  an  air  pump 
to  provide  oxygen  to  the  cathode,  whose  dynamics  are  denoted 
by  the  block  Ga  that  describes  the  relationship  between  the  air¬ 
flow  rate  N a  and  air  pump  voltage  VA.  The  coupled  blocks  G;, 
i=  1,  2,  3,  4,  describe  the  relationships  between  the  cell  current 
Ic  and  cell  voltage  Vc,  and  the  air-flow  rate  Aa  and  hydrogen 
flow  rate  Ah,  where  R  denotes  the  internal  resistance  of  the  fuel 
cell.  It  is  noticed  that  the  temperature  and  humidity  of  air  and 
hydrogen  are  not  managed  in  this  research,  but  are  regarded  as 
two  dependent  states  of  the  fuel  cell.  Assuming  that  the  fuel  cell 
is  operated  at  a  certain  operating  point,  the  six  blocks  can  be 
linearized  as  transfer  functions  of  finite  orders  and  time-varying 
coefficients.  These  coefficients  are  implicit  functions  of  state 
variables,  such  as  temperature,  membrane  hydration,  reactant 
pressure,  gas  concentration,  etc.,  which  directly  affect  cell  volt¬ 
age,  current  density,  and  the  overall  performance  of  the  fuel  cell. 
When  the  air  pump  voltage  and  the  hydrogen  valve  voltage  are 
regarded  as  input  variables,  and  the  cell  voltage  Vc  and  current 
Ic  are  chosen  as  system  outputs,  the  fuel  cell  system  model  can 
be  expressed  as 

Vc  =  Ga2yA +  GS4 VS -7?/c,  (1) 

/c  =  GaiVA  +  Gs3Vs,  (2) 

where  Gai=GaGi,  Ga2  =  GaG2,  Gs3  =  GsGs,  Gs4  =  GsG4. 
Since  Vc  and  Ic  are  dependent,  (1)  can  be  reduced  to 

Vc  =  (Ga2  -  RG^)Va  +  (GS4  -  RGM.  (3) 

It  is  noticed  that  the  cell  voltage  and  cell  current,  since  depen¬ 
dent,  are  not  to  be  regulated  concurrently  with  arbitrary  set 
points.  Either  configuration  (2)  or  (3),  namely  a  multi-input 
single-output  model,  is  employed  as  the  basis  of  system  iden¬ 
tification,  as  well  as  an  adjustable  model  for  designing  the  fuel 
cell  controller. 


Table  1 

Operating  conditions  for  system  identification 
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Load  8.5  V 


621 


Temperature  (°C) 

17 

Pressure  (atm) 

1 

Humidity  (%) 

74 

Fuel  cell  temperature  (°C) 

35-40 

Air  pump  voltage  (V) 

5  +  2 

Anode  reactant 

H2,  1.2  LPM,  7  psig 

Cathode  reactant 

Natural  air 

Current  density  (A  cm-2) 

0.2 

Sampling  frequency  (Hz) 

100 

2.3.  Off-line  system  identification 


The  purpose  of  off-line  system  identification  is  to  obtain  the 
transfer  functions  between  system  inputs  and  outputs  in  (1)  and 
(2)  at  certain  operating  points.  The  information  provided  from 
transfer  functions,  such  as  system  order,  parameters,  and  time 
delay,  must  be  used  for  the  adaptive  controller  design.  Although 
many  approaches  are  available,  non-parametric  system  identifi¬ 
cation  is  one  of  the  most  effective  to  investigate  the  frequency 
response  of  the  fuel  cell  stack. 

The  identification  of  the  transfer  function  Gai  =  IJV a  is 
demonstrated  between  the  cell  current  and  air  pump  voltage 
under  certain  operating  conditions,  as  described  in  Table  1.  The 
hydrogen  is  supplied  with  a  nominal  rate  of  1 .2  LPM  at  7  psig  by 
a  solenoid  valve,  which  is  always  on  during  the  process,  but  the 
output  valve  opens  to  purge  chemical  product  for  2  s  every  2  min. 
Even  though  the  purge  valve  is  closed,  the  hydrogen  is  contin¬ 
ually  supplied  and  consumed  to  convert  to  electricity,  while  the 
water  purging  can  be  regarded  as  the  external  disturbance  that 
helps  excite  system  dynamics  during  identification.  On  the  other 
hand,  this  disturbance  can  be  easily  controlled  by  the  proposed 
adaptive  controller. 

The  cell  voltage  is  set  fixed  by  an  electronic  load  meter  at 
8.5  V  for  three  or  more  identifications.  The  input  voltage  of  the 
air  pump  is  a  sinusoidal  function  Va  =  5  +  2  cos  (cot),  whose  fre¬ 
quency  sweeps  from  0.05  to  10  Hz.  Then,  the  dynamic  signal 
analyzer  (DSA)  HP35670A  takes  N  samples  of  cell  current  Ic 
and  calculates  the  magnitude  and  phase  of  frequency  response 
Gai(j&0  as  follows. 

|Gal(j«)l  =  \]  Jf  +  /c>  (4) 


<+(j<w)  =  arg  Gai(j®)  =  —  arctan 


(5) 


where 

1  n  j  N 

fs  =  —  y^/c(0  sin(<£>0  and  fc  =  —  F/c(0  cos(wf).  (6) 

t=  1  t=  1 

The  corresponding  frequency  responses  of  Gai  are  shown  in 
Fig.  4,  where  three  test  results  do  not  have  significant  differences 
over  the  frequency  range,  and  the  resulting  estimation  of  transfer 
function  is  then  expressed  in  different  orders  as 


io'1  io°  io1 

frequency  (Hz) 

Fig.  4.  Identified  frequency  responses  of  Ga i  at  a  load  of  8.5  V  for  three  tests. 


First-order  approximation: 

Ic(s)  0.01444s  +  0.08766 

Ga1  (s)  =  -  =  - 

w  Va(s)  s  + 0.229 

Second  order  approximation: 


Gal  (s)  = 


0.0004543s2  +  0.1307s  +  0.05606 
s2  +  1.011s  +  0.1392 


Third  order  approximation: 


(7) 


(8) 


Gai(s)  = 


0.002546s3  +  0.1356s2  +  0.06633s  +  0.02382 
s3  +  1.058s2  +  0.4321s  +  0.06131 


(9) 


Similar  results  can  be  also  obtained  for  the  transfer  func¬ 
tions  Ga 2,  Gs3,  and  Gs4,  and  bring  the  conclusion  that  a 
second-order  transfer  function  may  be  sufficient  to  describe 
the  dominant  dynamics  of  a  PEM  fuel  cell  stack  [18,19]. 
To  facilitate  the  following  adaptive  controller  design,  a  set 
of  second-order  continuous-time  transfer  functions  and  their 
discrete-time  approximations  of  0.01  s  sampling  time  at  a  certain 
operating  point  are  illustrated  in  Table  2. 


Table  2 


Illustration  of  identified  second-order  transfer  function  pairs  in  s  and  z  domains 


Cal  =Ic/Va,  Vc  =  8  V 

Ca2  =  Vc/V\,  7C  =  3  A 

0.9909s  +  15.89 

0.3269s  +  2.016 

s2  +  6.401s  +  9.732 

0.01037z  -  0.008835 

s2  +  3.921s  +  2.058 

0.003305z  -  0.003107 

z2  -  1.937z  + 0.938 

z2-  1.961z  +  0.9616 

GS3=Ws,  Vc  =  8V 

Gs4  =  Vc/Vs,/c=3A 

0.2477s  +  3.973 

0.08173s +  0.5041 

s2  +  6.401s  +  9.732 

0.002594z  -  0.002209 

s2  +  3.921s  +  2.058 

0.0008263z  -  0.0007769 

z2  -  1.937z  + 0.938 

z2-  1.961z  +  0.9616 

622 


Y.-P.  Yang  et  al.  /  Journal  of  Power  Sources  179  (2008)  618-630 


3.  Theory  and  calculation 

3.1.  On-line  parameter  estimation  and  adaptive  control 

PEM  fuel  cells  have  a  number  of  attributes  to  make  them 
promising  candidates  for  use  in  domestic  appliances  and  trans¬ 
portation  applications  [25].  The  power  demand  of  automobiles 
depends  on  the  driving  conditions,  and  the  provision  of  elec¬ 
tricity  from  a  fuel  cell  varies  due  to  its  system  states,  such  as 
reactant  flow  rates,  cell  temperature,  humidity,  and  pressure.  The 
relationship  between  fuel  cell  output  voltage  and  current  is  usu¬ 
ally  described  by  polarization  curves,  which  also  vary  with  cell 
temperature,  pressure,  humidity,  and  some  other  states.  These 
time-varying  and  state-dependent  properties  are  hardly  depicted 
by  simple  mathematical  models,  but  the  system  identification 
proves  that  it  is  possible  to  approximate  fuel  cell  dynamics  by 
transfer  functions  at  a  certain  operating  point.  These  transfer 
functions  can  be  transformed  to  discrete-time  parameter  estima¬ 
tion  models  in  the  adaptive  control  loop  of  a  fuel  cell  system 
subject  to  external  disturbances  as  well  as  system  parameter 
changes. 

In  practice,  the  fuel  cell  current  can  be  regarded  as  an  external 
disturbance  to  the  system,  which  is  requested  by  the  demand 
for  power  by  external  loads.  Therefore,  the  approximated  block 
diagram  of  a  fuel  cell  in  Fig.  3  can  further  simplified  as  a  two- 
input  and  one-output  structure. 

3.2.  ARM  A  model  and  on-line  parameter  estimation 

As  mentioned  above,  cell  voltage  and  cell  current  are  depen¬ 
dent,  and  thus  are  not  to  be  regulated  concurrently  with  arbitrary 
set  points.  Either  (2)  or  (3)  is  used  for  on-line  parameter  esti¬ 
mation,  and  is  usually  represented  as  a  discrete-time  equation 
of  difference  operator,  leading  to  the  following  auto-regressive 
moving-average  (ARMA)  model  [26] : 

A(q ,  k)y(k)  =  B(q,  k)u(k)  +  w(k),  (10) 

A(q,  k)  =  ao(k)  +  a\(k)q~l  +  •  •  •  +  ar(k)q~r ,  (11) 

B(q,  k)  =  Bx(k)q-1  +  B2(k)q~2  +  •  •  •  +  Br(k)q~r ,  (12) 

where  y  =  Ic  or  Vc  denotes  the  output  variable,  u  =  [u\  u2]T  =  [Va 
Vs]T  represents  the  2-tuple  input  vector,  w  is  composed  of 
higher-order  dynamics  and  disturbances,  and  is  known  as  an 
estimation  error  vector.  The  coefficients  af  s  are  scalars  while 
Bi  s  are  2- tuple  row  vectors  for  i=  1,  2,  . . .,  r,  the  integer  k 
denotes  the  time  instant  of  sampling,  and  r  is  the  system  order. 
Without  loss  of  generality,  ao  =  1 .  In  a  succinct  form,  the  term 
A(q,k)y(k)  is  the  auto-regression  part,  where  q~l  is  a  backward 
shift  operator,  and  q~hy(k)  represents  the  past  output  y(k  —  h ); 
h-  1 ,  2,  ...  r;  B(q,k)u(k)  denotes  the  moving  average  of  past 
inputs,  where  q~hu(k )  =  u(k  —  h). 

The  approach  of  the  recursive  least  square  algorithm  for  sys¬ 
tem  identification  is  to  predict  the  system  output  according  to 
the  past  information  of  input  and  output  measurements,  as  well 
as  the  updated  set  of  system  parameters.  From  (10)  to  (12),  the 


estimation  of  y(k)  is  expressed  as 

%k)  =  4>J(k  -  1)  9(k  -  1),  (13) 

where  <fiT(k—  1)  =  [— y(k—  1)  ...  —  y(k  —  r)  u\(k—l)  ... 
u\(k  —  f)  u2(k  —  1)  . . .  u2{k  —  r)]  and  6{k  —  1)  is  the  estimated 
parameter  vector  of  a't s  and  B\ s.  The  estimation  error  w(k)  is 
defined  as  y(k)  —  y(k). 

At  step  k ,  when  the  new  output  is  measured,  the  updated  set 
of  parameters  are  calculated  by  [26] 

6(k)  =  6(k  -  1)  +  K(k)[y(k)  -  4>T(k  -  1  )6(k  -  1)],  (14) 

K(k )  =  P(k  -  1  )<t>{k  -  1)[A  +  ct7(k  -  1  )P(k  -  1  )<t>(k  -  1)]~\ 

(15) 

n/(1  [I-K(k)<t?(k-\)]P(k-\) 

P(k)  = - - - ,  (16) 

where  the  estimation  gain  K(k)  brings  the  relative  information 
of  new  measurements  to  update  the  parameter  estimation,  and 
the  covariance  matrix  P(k)  characterizes  the  difference  between 
the  estimated  parameters  and  their  true  values.  Initially,  P( 0) 
is  chosen  as  a  large  positive  number,  for  example  103/,  where 
I  is  an  identity  matrix,  to  render  the  inaccurate  initial  guess 
§(0)  quickly  negligible.  The  coefficient  0  <  k  <  1  is  called  the 
forgetting  factor,  which  places  more  importance  on  the  new 
information  for  updating  system  parameters,  and  less  attention 
on  the  past  information. 

3.3.  Model  reference  adaptive  control  (MRAC) 

Model  reference  adaptive  control  is  one  of  the  promising 
strategies  for  a  complex  fuel  cell  system.  No  matter  how  com¬ 
plicated  the  system  is,  its  desired  output  can  be  designated  to 
follow  the  output  of  a  reference  model  with  specified  dynam¬ 
ics.  The  model  reference  adaptive  control  strategy  includes  the 
specification  of  a  reference  model  with  the  desired  dynamics, 
on-line  parameter  estimation,  and  calculation  of  control  gains. 
This  adaptive  control  scenario  is  illustrated  in  Fig.  5. 

From  the  previous  result  of  system  identification,  the  approx¬ 
imated  blocks  of  the  fuel  cell  system  of  interest  can  be  well 
described  as  second-order  subsystems.  Without  loss  of  gener- 


Fig.  5.  Model  reference  adaptive  control  scenario  of  fuel  cell. 
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ality,  the  following  control  law  is  derived  for  an  SISO  system 
whose  model  is  identified  on-line  and  is  denoted  by 

A{q)y{k)  =  B{q)u(k),  (17) 

where  A(q)  =  q2  +  a\q  +  a2 ,  and  B(q)  =  b\q  +  b 2.  For  a  stan¬ 
dard  feedback  control  law 


R(q)u(k)  =  T(q)uc(k)  +  S(q)y(k), 


(18) 


where  uc  is  a  reference  signal.  The  closed-loop  relationship 
between  uc(k)  and  y(k)  can  be  written  as 


y(k)  = 


BT 

AR  +  BS 


uc(k ), 


(19) 


in  which  the  index  q  is  omitted  for  simplicity.  The  desired  closed- 
loop  model  is  chosen  as 


yi 


6- 


y 


yj 


Fig.  6.  Multivariable  adaptive  control  block  diagram. 


and  air  pump,  which  accounts  for  the  nominal  ratio  between  the 
flow  rates  of  hydrogen  and  oxygen.  The  determination  of  control 
functions  R,  T,  and  S  still  applies  through  (23)-(27)  by  simply 
replacing  B(q)  by  B\(q)  +  KB2(q). 


—  BmUc  (k), 


(20) 


where  Am(q)  =  q2  +  am  1  q  +  am 2  and  Bm(q)  =  bm\  q  +  bm 2,  then  a 
perfect  model-following  condition  becomes 


BT  =  B^ 
AR  T  BS  Am 


(21) 


Under  the  conditions  of  compatibility  and  causality,  the 
control  polynomials  T  and  S  are  obtained  by  the  Bezout  or 
Diophantine  equation: 


AR  T  BS  —  AmA0, 


(22) 


where  R(q)  =  q  +  r\ ,  S(q)  =  soq  +  sq ,  and  A0(q)  =  q  +  a0  is  known 
as  an  observer  polynomial,  which  is  always  cancelled  in  the 
transfer  function  between  the  output  and  reference  signal.  When 
the  unit  steady- state  gain  of  Bm/Am  is  designed  and  Bm  =  /3B  is 
selected,  the  control  polynomial  T  becomes  simply  /3A0.  In  sum¬ 
mary,  the  parameters  of  those  control  polynomials  are  calculated 
after  system  parameters  a\,  a2 ,  b\,  and  b2  are  estimated: 


n= 


(A2  *fin2  ^o^ml)^1^2  (^o  T  ^ml  Cli^b2 


(23) 


4.  Result  and  discussion 

From  the  perspective  of  the  certainty  equivalence  principle  in 
the  system  identification  and  adaptive  control  theories,  a  linear 
model  of  reduced  order  with  time- varying  parameters  is  usually 
used  to  generate  the  control  law  as  if  it  were  a  true  system  no 
matter  whether  it  is  linear  or  nonlinear  [26] .  Five  adaptive  con¬ 
trol  scenarios  are  investigated:  the  regulation  of  cell  voltage  by 
the  flow  rate  of  air  or  hydrogen,  the  regulation  of  cell  current 
by  the  flow  rate  of  air  or  hydrogen,  and  the  regulation  of  cell 
voltage  by  the  flow  rates  of  air  and  hydrogen.  These  scenar¬ 
ios  provide  a  systematic  way  of  demonstrating  the  effectiveness 
of  the  proposed  controller  for  various  situations  of  a  change  in 
the  external  load  with  which  the  change  of  internal  impedance 
and  the  corresponding  states  accompany.  The  first  four  scenarios 
are  SISO  cases,  while  the  last  one  pertains  to  the  multivariable 
configuration. 

The  experimental  setup  of  the  fuel  cell  stack  and  its  periph¬ 
eral  devices  is  shown  in  Fig.  7.  On-line  parameter  estimation 
and  model  reference  adaptive  control  law  are  performed  on  a 
laptop  computer  using  Matlab.  The  fuel  cell  power  is  dissi¬ 
pated  on  an  electronic  load  meter  by  which  either  cell  current  or 


so  = 


b2(a0am\  -CL2-  am\a\  +  am2  -  a\a0)  +  b\{am\a2  -  a\a2  -  a0am2  -  a0a2) 


si  = 


b2(a\a2  -  am\a2  +  a0am2  -  a2a0)  +  b\(am2a2  -  a2  -  a0am2a\  -  a0a2am\) 


(24) 

(25) 


ot  =  b  2  —  a\b\b>2  +  a2b\. 


(26) 


bm\  +  bm2 
b\  +  b>2 


(27) 


In  the  following  experiments,  am i=0,  am2  =  0.01,  and 
a0  =  —  0.1,  which  correspond  to  the  desired  closed-loop  poles 
assigned  at  —0.1  and  itO.ly'  within  the  unit  circle  of  z-plane. 
The  above  algorithm  for  MRAC  can  be  extended  to  the  case 
of  two  inputs,  air  and  hydrogen  flow  rates,  and  one  output,  cell 
voltage,  as  shown  in  Fig.  6.  The  parameter  k  is  defined  as  the 
nominal  ratio  between  the  average  voltages  of  solenoid  valve 


cell  voltage  remains  constant  while  the  other  is  regulated.  Their 
dependence,  however,  is  characterized  by  polarization  curves  as 
functions  of  pressure,  humidity,  and  temperature.  The  air-flow 
rate  is  controlled  by  an  air  pump,  whose  property  is  specified 
by  experiments  with  the  relationship  of  flow  rate  and  driving 
voltage,  as  shown  in  Fig.  8.  The  hydrogen  flow  rate  is  kept  con¬ 
stant  at  1.2  LPM  at  7  psig  in  the  SISO  cases,  but  is  adaptively 
controlled  by  the  solenoid  valve  at  on-off  sequences.  Both  the 
air  and  hydrogen  are  not  humidified,  and  the  cell  temperature  is 
monitored  and  sequentially  controlled  around  35-40  °C. 
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Fig.  7.  (a)  Fuel  cell  control  block  diagram  and  (b)  experimental  setup. 


Load  4A  ->  3A  ->  2A 
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Fig.  8.  Air-flow  rate  vs.  air  pump  voltage. 


Fig.  9.  Cell  voltage  control  with  fixed  system  parameters. 


Y.-P.  Yang  et  al.  /  Journal  of  Power  Sources  179  (2008)  618-630 


625 


(b)  Load  8V 
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Fig.  10.  (a)  Adaptive  control  of  cell  current  (2  A)  by  the  air  pump  voltage  (upper 
curve)  under  the  cell  voltage  set  at  8  V,  and  (b)  the  corresponding  parameter 
estimation. 


Fig.  11.  (a)  Adaptive  control  of  cell  current  (solid  curve)  at  3  A  for  the  cell 
voltage  (dashed  line)  of  9  V,  8.5  V,  and  8  V,  (b)  the  input  air  pump  voltage,  and 
(c)  the  corresponding  parameter  estimation. 


It  is  noteworthy  that  the  fuel  cell  has  a  built-in  water  and 
thermal  management  subsystem;  thus,  all  the  following  experi¬ 
ments  focus  on  the  justification  and  effectiveness  of  the  proposed 
adaptive  control  strategy.  In  addition,  the  fuel  cell  used  in  this 
research  was  not  in  its  best  condition  due  to  the  damage  and 
aging  of  membranes,  so  that  the  fuel  stack  could  only  oper¬ 
ate  in  a  relatively  conservative  range  of  voltage,  current,  and 
power. 

4.1.  Control  of  cell  voltage  with  fixed  parameters 

Before  the  adaptive  control  performance  is  investigated,  it 
is  interesting  to  observe  the  performance  of  cell  voltage  that 
is  adjusted  by  air  pump  voltage  according  to  the  control  law 
with  fixed  control  parameters  in  (23)— (27).  These  parameters  are 
determined  by  a  fixed  set  of  system  parameters  a\  =  —1.1059, 
b\  =  8.7898,  <22  =  0.3567,  and  Z?2  =  0.2405,  while  the  hydrogen  is 
supplied  at  a  constant  flow  rate  of  1 .2  LPM  at  7  psig.  As  shown  in 


Fig.  12.  (a)  Adaptive  control  of  cell  voltage  (solid  curve)  at  8.5  V  under  the 
load  demand  (dashed  line)  of  4  A,  3  A,  and  2  A,  (b)  the  corresponding  air  pump 
voltage,  and  (c)  parameter  estimation. 
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Fig.  9,  the  cell  voltage  does  not  approach  the  desired  command 
of  9  V  subject  to  the  variation  of  the  current  load  from  4  A  and 
3  A  to  2  A  for  200  s  each.  It  is  concluded  that  the  conventional 
control  strategy  derived  from  fixed  system  parameters  seems  to 
lack  robustness  to  external  disturbances  as  well  as  to  the  variation 
of  system  parameters,  and  the  large  variation  in  cell  voltage  will 
cause  the  DC/DC  converter  to  regulate  the  output  voltage  with 
lower  efficiency. 

4.2.  Regulation  of  cell  current  by  air-flow  rate 

First,  adaptive  control  of  the  fuel  cell  current  by  the  air 
pump  voltage  corresponding  to  the  air  pump  flow  rate  is  demon¬ 
strated.  When  the  output  voltage  is  fixed  by  the  load  meter  and 
the  hydrogen  flow  is  kept  steady  by  the  solenoid  valve,  the 
cell  current  is  adaptively  adjusted  by  the  air  pump  voltage  on 
the  demand  for  power  consumption  subject  to  the  plant  vari¬ 
ation.  This  SISO  process  is  simply  described  by  the  transfer 
function  Gai  which  is  transformed  into  a  discrete-time  model 
and  is  easily  expressed  as  an  ARMA  model  with  time- varying 
parameters. 

For  a  2  A  current  demand,  the  load  voltage  is  set  at  8  V,  and 
the  hydrogen  is  supplied  continuously  with  1.2  LPM  at  7psig. 
Notice  that  hydrogen  is  oversupplied  under  this  situation.  The 
time  responses  of  the  air  pump  voltage  input  and  cell  current  out¬ 
put  are  shown  in  Fig.  10(a),  while  the  corresponding  variation 
of  parameters  estimated  on-line  is  depicted  in  Fig.  10(b).  It  is 
found  that  the  system  parameters  did  not  vary  much,  and  the  cell 


current  was  easily  regulated  within  2  zb  0.02  A  rms  by  adjusting 
the  air-flow  rate.  When  the  cell  voltage  shifts  from  9  V  through 
8.5  V  to  8  V,  at  200  and  400  s,  the  current  is  regulated  at  3  A,  as 
shown  in  Fig.  11(a),  by  adaptively  adjusting  the  air  pump  volt¬ 
age  to  change  the  air-flow  rate  in  Fig.  1 1(b).  The  corresponding 
parameter  estimation  is  displayed  in  Fig.  1 1(c).  Various  settings 
of  cell  voltage  mean  different  operating  points,  which  result  in 
different  parameters  in  the  second-order  transfer  function.  At 
the  junctions  of  voltage  change,  the  cell  current  experienced 
abrupt  increase,  but  was  immediately  adapted  to  its  desired 
value  within  ±0.02  A  rms.  This  experiment  shows  that  the  pro¬ 
posed  adaptive  control  is  robust  to  the  plant  variation  and  power 
demand. 

4.3.  Regulation  of  cell  voltage  by  air-flow  rate 

In  many  applications,  the  power  source  has  to  be  kept  at  a 
constant  voltage,  while  the  current  is  drawn  based  on  the  demand 
of  the  load.  This  happens  quite  often  for  fuel  cells  used  on  electric 
vehicles.  If  hydrogen  is  constantly  oversupplied  with  1 .2  LPM  at 
7  psig,  and  the  cell  current  is  drawn  from  4  A  through  3  A  to  2  A, 
the  voltage  will  be  well  regulated  at  8.5  V  by  the  air-flow  rate, 
as  shown  in  Fig.  12(a).  The  corresponding  time  history  of  air 
pump  voltage  is  presented  in  Fig.  12(b).  It  is  not  surprising  that 
the  system  parameters  in  Fig.  12(c)  vary  but  are  soon  identified 
at  different  operating  points,  while  the  control  law  is  adaptively 
changed  to  keep  the  cell  voltage  within  a  variation  of  ±0.1  V  of 
the  desired  value. 


Load  8.5V 
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Fig.  13.  (a)  Adaptive  control  of  cell  current  by  (b)  the  voltage  signal  to  the  solenoid  valve  with  (c)  its  corresponding  hydrogen  flow  rate,  while  (d)  system  parameters 
are  estimated. 
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4.4.  Regulation  of  cell  current  by  hydrogen  flow  rate 

It  is  also  interesting  to  investigate  how  well  the  cell  current 
is  controlled  by  the  hydrogen  flow.  When  the  air  pump  voltage 
is  set  at  5  V  with  about  6. 6-6. 8  LPM  air-flow,  the  cell  voltage 
is  fixed  at  8.5  V,  and  the  desired  current  is  4  A,  the  output  cur¬ 
rent  approaches  4.08  A  with  a  minor  error  of  2%,  as  shown  in 
Fig.  1 3(a).  Although  the  cell  current  presents  minor  oscillation  of 
±0.08  A,  the  hydrogen  flow  rate  is  measured  about  0.123  LPM, 
which  is  only  about  10%  consumption  compared  to  the  exces¬ 
sive  usage  of  hydrogen  at  1 .2  LPM.  During  the  adaptive  control 
process,  a  PWM  circuit  is  designed  to  control  the  solenoid  valve 
at  100  Hz.  The  voltage  signal  controlling  the  on-off  actuation  of 
solenoid  valve  is  shown  in  Fig.  13(b),  and  is  used  to  calculate 
the  average  hydrogen  flow  rate  over  each  second,  as  shown  in 
Fig.  13(c).  The  corresponding  on-line  parameter  estimates  are 
displayed  in  Fig.  13(d)  for  the  linearly  approximated  system 
model. 

With  the  same  setting  of  air  pump  voltage,  the  cell  voltage  is 
fixed  at  8.5  V  for  the  first  200  s,  9  V  for  the  next  200  s,  and  9.5  V 
for  the  last  200  s,  while  the  desired  cell  current  is  specified  at 
5  A,  4  A  and  3  A.  As  shown  in  Fig.  14(a),  the  cell  current  is  con¬ 
trolled  successfully  by  adjusting  the  flow  rate  of  hydrogen  with  a 
solenoid  valve  as  depicted  in  Fig.  14(b),  and  the  corresponding 
parameter  estimation  is  illustrated  in  Fig.  14(c).  As  expected, 
77.3%  of  hydrogen  is  saved  during  the  first  200  s,  84%  during 
the  next  200  s,  and  85%  for  the  last  200  s,  compared  with  the 
constant  hydrogen  flow  of  1.2  LPM.  It  is  also  interesting  to  find 
that  the  steady-state  error  of  current  response  is  within  ±0.35% 
with  minor  oscillations. 

4.5.  Regulation  of  cell  voltage  by  hydrogen  flow  rate 

In  this  experiment,  the  air  pump  voltage  is  fixed  at  5  V  to 
which  the  corresponding  flow  rate  is  around  6.6-6. 8  LPM,  while 
the  hydrogen  pressure  is  kept  at  7  psig.  The  operating  state  of 
the  solenoid  valve  is  either  fully  on  at  12  V  or  fully  off  at  0  V. 
The  flow  rate  of  hydrogen  can  be  controlled  by  adjusting  the 
duty  cycle,  or  the  open  period  of  the  valve,  which  is  determined 
adaptively  by  the  control  law  and  a  carrier  signal  to  chop  into 
a  pulse  width  modulation  (PWM)  signal.  Fig.  15(a)  shows  the 
controlled  cell  voltages  from  9.5  V  through  9  V  to  8.5  V  on  the 
demand  of  a  3  A  load,  controlled  by  the  hydrogen  flow  rate  in 
Fig.  15(b),  and  the  corresponding  time- varying  parameters  are 
estimated  in  Fig.  15(c). 

It  is  found  that,  at  different  power  demands,  the  cell  voltage 
controlled  by  the  hydrogen  flow  presents  larger  steady- state  error 
and  larger  amplitude  of  oscillation  that  are  controlled  by  the 
air-flow.  In  the  first  200  s,  the  average  voltage  is  9.54  V  with 
0.43%  error,  and  the  oscillatory  amplitude  is  within  ±0.125  V; 
for  the  next  200  s,  the  average  voltage  is  9. 1 8  V  with  2%  error  and 
±0.225  V  oscillation;  while  the  average  voltage  is  8.78  V  with 
3.26%  error  and  ±0.375  V  oscillation  for  the  last  200  s.  Higher 
oscillation  of  cell  voltage  occurs  at  a  lower  power  output  because 
of  the  pulsation  and  overcharge  of  hydrogen  flow  due  to  the 
on-off  action  of  the  solenoid  valve.  This  can  be  improved  by  the 
use  of  a  proportional  solenoid  valve  whose  opening  can  be  finely 
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Fig.  14.  (a)  Adaptive  control  of  cell  current  by  (b)  hydrogen  flow  rate  while  (c) 
system  parameters  are  estimated. 

adjusted.  However,  the  amount  of  hydrogen  consumed  is  much 
less  than  the  continuous  feeding  of  1.2  LPM.  From  Fig.  15(b), 
the  average  flow  rate  of  hydrogen  is  0.17  LPM  during  the  first 
200  s,  and  86.2%  hydrogen  is  saved;  the  average  flow  rate  of 
hydrogen  is  0.1 6  LPM  during  the  next  200  s  with  a  saving  of 
86.9%;  while  88.6%  hydrogen  is  saved  with  the  flow  rate  at 
0.14  LPM  during  the  last  200  s. 

In  the  case  when  the  cell  voltage  is  controlled  at  9  V  on  var¬ 
ious  demands  of  load  from  4  A  through  3  A  to  2  A,  as  shown 
in  Fig.  16(a),  the  response  of  the  cell  voltage  also  shows  an 
oscillatory  behavior,  especially  at  lower  power  due  to  the  reason 
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time  (s) 


Fig.  15.  (a)  Adaptive  control  of  cell  voltage  from  9.5  V  through  9  V  to  8.5  V 
at  3  A  load,  (b)  the  required  flow  rate  of  hydrogen,  and  (c)  the  corresponding 
parameter  estimation. 

mentioned  in  the  previous  case.  The  average  voltage  is  9.08  V 
with  0.92%  error  for  the  first  200  s,  and  the  average  voltage 
is  9.11  V  with  1.26%  error  for  the  next  200  s,  while  the  aver¬ 
age  voltage  is  9.34  V  with  3.74%  error  for  the  last  200  s.  The 
corresponding  input  of  the  hydrogen  flow  rate  and  the  parame¬ 
ter  variation  are  shown,  respectively,  in  Fig.  16(b)  and  (c).  The 
hydrogen  consumption  is  saved  for  84.8%  for  4  A  load,  90.3% 
for  3  A  load,  and  92.9%  for  2  A  load,  in  comparison  with  the 
constant  flow  rate  of  1.2  LPM. 

4.6.  Regulation  of  cell  voltage  both  by  air  and  hydrogen 
flow  rates 

In  the  previous  cases,  the  flow  rate  of  either  air  or  hydro¬ 
gen  is  able  to  control  the  cell  voltage  and  current  adaptively 
based  on  the  SISO  model.  Would  the  air  pump  and  solenoid 
valve  work  together  to  regulate  the  cell  voltage  due  to  the 
demand  of  the  load  current?  Theoretically,  the  stoichiometry 


(c)  Load  4A  ->  3A  ->2A 


Fig.  16.  Adaptive  control  of  cell  voltage  at  load  from  4  A  through  3  A  to  2  A, 
(b)  the  required  flow  rate  of  hydrogen,  and  (c)  the  corresponding  parameter 
estimation. 

of  hydrogen  and  oxygen  in  the  chemical  reaction  of  PEM  fuel 
cell  is  2:1.  If  the  cathode  side  is  supplied  with  air,  the  stoichio¬ 
metric  ratio  between  hydrogen  and  air  becomes  approximately 
2:5.  In  practice,  an  excessive  amount  of  oxygen  is  requested 
for  the  fuel  cell  to  operate  stably  [15].  From  the  above  exper¬ 
iments,  the  parameter  k  that  denotes  the  ratio  between  the 
average  voltages  of  the  solenoid  valve  and  air  pump  was  found 
at  a  nominal  number  of  0.04  as  the  fuel  cell  operates  in  a 
steady-state.  This  nominal  number  will  be  applied  for  the  mul¬ 
tivariable  adaptive  control  described  in  Fig.  6  in  the  following 
experiments. 

The  pressure  of  hydrogen  from  the  storage  still  maintains  at 
7  psig  before  entering  the  fuel  cell  through  the  solenoid  valve. 
Both  the  air  pump  and  solenoid  valve  supply  oxygen  and  hydro¬ 
gen  with  a  proportional  ratio  of  k  to  regulate  the  cell  voltage  at 
9.5  V  while  the  load  current  demand  varies  from  4  A  through 
3  A  to  2  A  at  200  and  400  s,  as  shown  in  Fig.  17(a).  The  adaptive 
control  law  of  voltages  of  the  air  pump  and  solenoid  valve  is 
shown  in  Fig.  17(b);  the  corresponding  hydrogen  flow  rate  is 
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Fig.  17.  (a)  Adaptive  control  of  cell  voltage  (solid  curve)  under  the  current 
command  (dashed  line),  (b)  the  control  voltage  of  air  pump  and  solenoid  valve, 
and  (c)  the  corresponding  parameters. 


displayed  in  Fig.  17(c),  and  the  estimation  curves  of  the  system 
parameters  are  illustrated  in  Fig.  17(d). 

Compared  to  the  continuous  hydrogen  supply  at  1.2  LPM, 
the  average  flow  rate  of  0.13  LPM  during  the  first  200  s  for 
the  4  A  load  accounts  for  a  reduction  of  89.9%  hydrogen  con¬ 
sumption.  During  the  next  200  s  for  the  3  A  load,  the  hydrogen 
consumption  is  reduced  to  90.6%  at  the  average  flow  rate  of 
0.1 14  LPM,  and  during  the  last  200  s  for  the  2  A  load,  94% 
hydrogen  is  saved  at  the  average  flow  rate  of  0.087  LPM.  At 
the  transition  of  load  variation,  the  cell  voltage  experiences  an 
apparent  oscillation,  but  settles  quickly  within  a  few  seconds 
due  to  the  robustness  of  adaptive  control.  Due  to  the  regulation 
of  cell  voltage  by  the  hydrogen  flow  rate,  even  though  much 


Load  3A 


50  100  150  200  250  300 


time  (s) 

Fig.  18.  Cell  voltage  under  (a)  an  excessive  supply  (1.2  LPM)  and  (b)  economic 
supply  (0.12  LPM)  of  hydrogen. 


hydrogen  is  saved,  the  amplitude  of  the  sustained  oscillation  of 
cell  voltage  is  much  greater  than  that  under  an  excessive  supply 
of  hydrogen.  Fig.  18(a)  shows  the  well-controlled  cell  voltage 
by  an  excessive  supply  of  hydrogen  of  1.2  LPM  at  the  3  A  load, 
but  Fig.  18(b)  shows  a  heavier  oscillation  of  cell  voltage  as  the 
hydrogen  is  adjusted  economically  at  0.12  LPM  by  the  rapid 
parameter  adjustment  of  the  adaptive  controller. 

In  summary,  from  the  above  experiments,  when  the  stoichio¬ 
metric  ratio  of  hydrogen  and  oxygen  is  larger  than  2,  i.e.,  the 
hydrogen  is  oversupplied,  both  the  cell  voltage  and  current  can 
be  controlled  by  the  air- flow  rate  with  only  minor  oscillations. 
When  the  oxygen  or  air  is  oversupplied  and  in  contrast  the  hydro¬ 
gen  is  supplied  in  an  economic  mode  by  the  adaptive  control  law, 
both  the  cell  current  and  voltage  regulated  by  the  hydrogen  flow 
rate  present  larger  oscillations.  Table  3  summarizes  the  qualita¬ 
tive  description  of  the  resulting  oscillation  of  cell  outputs  to  the 
under-  or  over-supply  of  hydrogen.  However,  such  oscillations 
can  be  suppressed  by  a  proper  choice  of  DC/DC  converter.  In 
short,  the  proposed  model  reference  adaptive  strategy  results  in  a 
robust  performance  to  the  parameter  variation  at  different  power 
demands. 


Table  3 


Qualitative  analysis  of  cell  outputs 


Input  flow 

Stoichiometry  of  H2/O2 

Cell  output 

Oscillatory 

Air 

>2 

Current 

Minor 

Air 

>2 

Voltage 

Minor 

Hydrogen 

<2 

Current 

Medium 

Hydrogen 

<2 

Voltage 

Large 

Air  +  hydrogen 

<2 

Voltage 

Large 
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5.  Conclusion 

Indirect  model  reference  adaptive  control  has  been  success¬ 
fully  applied  to  a  low-power  PEM  fuel  cell.  For  use  in  control 
systems,  the  highly  nonlinear,  time- varying,  distributed  param¬ 
eter  fuel  cell  stack  of  complex  electrochemical  reactions  and 
physics  is  approximated  by  a  simple  two-input  two-output  model 
that  is  described  within  subsystem  blocks.  These  blocks  describe 
the  dynamic  relationship  between  cell  voltage  and  current  out¬ 
puts  and  the  hydrogen  and  air-flow  rate  inputs.  By  including  the 
dynamics  of  the  air  pump  and  solenoid  valve,  the  correspond¬ 
ing  transfer  functions  are  identified  off-line  at  certain  operating 
points  via  a  non-parametric  frequency  domain  method.  The 
model  is  subsequently  reduced  to  an  auto-regressive  moving- 
average  model  for  on-line  parameter  estimation  in  an  adaptive 
control  loop.  Experimental  results  show  that  the  adaptive  control 
is  effective  in  regulating  cell  current  and  voltage  by  adjusting 
air  and  hydrogen  flow  rates  in  response  to  fuel  cell  and  exter¬ 
nal  power  demand  perturbations.  On-off  control  of  the  hydrogen 
valve  is  shown  to  reduce  fuel  consumption,  which  increases  sys¬ 
tem  efficiency.  It  is  noteworthy  that  the  successful  regulation  of 
cell  voltage  by  manipulation  of  the  air  pump  and  hydrogen  valve 
may  relieve  DC/DC  converters  of  demanding  requirements  for 
regulating  fluctuating  cell  voltage.  Nevertheless,  the  model  ref¬ 
erence  adaptive  control  is  applied  under  restricted  conditions  of 
operation  in  this  research;  it  is  necessary  to  make  more  effort 
toward  its  robustness  and  generality  to  a  class  of  fuel  cell  sys¬ 
tems.  In  addition,  it  is  an  important  issue  to  include  the  water  and 
thermal  management  in  future  research  to  monitor  and  control 
the  humidity  and  temperature  of  reactants  and  the  membrane 
so  that  the  fuel  cell  can  operate  in  a  safe,  efficient,  and  robust 
way  with  the  proposed  control  strategy.  It  is  expected  that  the 
proposed  adaptive  control  can  be  easily  extended  to  more  com¬ 
plex  multivariable  control  problems  when  more  system  states 
are  to  be  regulated,  and  could  be  applied  to  stationary  and  trans¬ 
portation  systems,  as  well  as  low-  and  high-power  fuel  cell 
applications. 
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